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ABSTRACT 



We repeat and extend the analysis of Eriksen et al 2004 and Hansen et al 2004 
testing the isotropy of the Cosmic Microwave Background (CMB) fluctuations. 
We find that the hemispherical power asymmetry previously reported for the 
largest scales £ = 2 — 40 extend to much smaller scales. In fact, for the full 
multipole range ^ = 2 — 600, significantly more power is found in the hemisphere 
centered at {9 = 107° ± 10°, (f) = 226° ± 10°) in galactic co-latitude and longitude 
than in the opposite hemisphere consistent with the previously detected direction 
of asymmetry for £ = 2 — 40. We adopt a model selection test where the direction 
and amplitude of asymmetry as well as the multipole range are free parameters. 
A model with an asymmetric distribution of power for £ = 2 — 600 is found to 
be preferred over the isotropic model at the 0.4% significance level taking into 
account the additional parameters required to describe it. A similar direction of 
asymmetry is found independently in all six subranges of 100 multipoles between 
£ = 2 — 600 and none of our 9800 isotropic simulated maps show a similarly 
consistent direction of asymmetry over such a large multipole range. No known 
systematic effects or foregrounds are found to be able to explain the asymmetry. 

Subject headings: (cosmology:) cosmic microwave background — cosmology: 
observations — methods: data analysis — methods: statistical 



Introduction 



In the first public release of the WMAP data, we reported a signifi c ant asymme- 



try of the distribution of large scale power on the sky (jEriksen et al.l l2004al : iHansen et al. 



2004al ). This finding as well as other statistical anomalies were confirmed by several other 



authors using other methods (de Oliveira-Costa et al. 


2004; 


Teemar 


k et al. 2003: Park 


2004; 


Vielva et al.ll2004b: 


Eriksen et al.ll2004b: 


jand & Maeueiio 


2005a b: 


Larson & Wandelt 


2004; 


McEwen et al. 2005 


: Hansen et a 


2004b 


c ) . Some of these have been confirmed again in the 


3 and 5 year data (] 


iinshaw et al. 


20071: J 


5pereel et al. 


2006 


: Cruz et al. 2006: McEwen et al. 


2006: Eriksen et all 


2007: Pietrobon et a 


L 2008). In f 


Hansen et al. 


2004a) we found that the 



independent multipole ranges ^ = 2 — 19 and £ = 20 — 40 were both particularly asymmetric 
but with two different axes of asymmetry, the former being a galactic north-south asymmetry 
and the latter being an east-west asymmetry. Using the full range £ = 2 — AO we found the 
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highest significance of the asymmetry with the axis pointing in the direction of (100°, 237°) 
in galactic co-latitude and longitude (this is the convention we will use for all sky positions 
in this paper). 



The position of the non-Gaussian cold spot found in (IVielva et al.ll2004bl ) is found to 
be positioned close to the center of the hemisphere with large fluctuation power. Further- 
more, the CMB signal also demonstrated non -Gaussian statistical properties in the direction 



where the power spectrum amplitude was low (jPark II2004J : lEriksen et al.ll2004bl : iHansen et al. 



2004bl ). Our own attempts to provide an explanation to the asymmetry invoked a class of 
homogeneous models that include anis otropic expansion (shear) and global rotation (vor- 



ticity) - the Bianchi type Vllh mode ls( Jaffe et al 



physical explanation (|Jaffe et al.ll2006l : [Bridges et al 



2005a|jb|). Unfortunately, this fails as a 



20071 ) since the best-fit parameter space 



is then inconsistent with a wide range of other evidence. 



Recent theoretical developments (I Ackerman et al.ll2007l : lErickcek et al.ll2008l : lGordon et al. 



20051 ) have proposed new mechanisms for generating the imprint of a preferred direction in 

the CMB. Some of these models would lead to hemispherical asymmetry also fo r multipoles 

much larger than i = 40 (angular scales smaller than 5 degrees). A recent analysis( lGroeneboom fc Eriksen 



20081 ) of the Ackerman et al. model indicates a good fit to the WMAP data up to a multipole 
moment i = 400. However, this model is not able to address the issue of the hemispherical 
power asymmetry. 

In this paper we reinvestiga te the nature of the power asymmetry using the recent 5 
year data release from WMAP (IHinshaw et al.ll2008l ). and specifically assess whether the 
signal manifests itself on smaller angular scales £ >> 40. Instead of randomly searching 
for an asymmetry in arbitrary multipole ranges, we adopt a model selection procedure that 
searches for the best fit asymmetry direction and multipole interval. An asymmetric model of 
the CMB needs more free parameters (i.e. direction of the axis of asymmetry) than a simpler 
isotropic model. The model selection procedure tests whether an asymmetric model of the 
CMB is actually preferred by the data taking into account the introduction of additional free 
parameters. We introduce a set of general asymmetric models with three to five additional 
parameters. These models are not based on any physical theory but are rather general 
parametrization of hemispherical asymmetry. We use these models in order to investigate 
how the asymmetry is distributed in harmonic space as well as on the sphere. 

In section [2] we describe the data and masks used in the analysis. The methods used 
to assess the asymmetry are outlined in section [31 In section S] we show the results of the 
isotropy tests applied to the 5 year WMAP data and in section Owe conclude. 
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2. Data 



The analysis in this paper was performed using the 5 year release of the WMAP data 
(publicly available at the Lambda web sit^) as well as a large ensemble of simulated maps of 
each channel Q (41GHz), V (61GHz) and W (94GHz) (the result for each channel is obtained 
by taking the mean of all DAs for each channel). From all maps, we have subtracted the best 
fit mono- and dipole. Some tests are also performed on maps from single DAs and single 
years of observation. 

A series of galactic masks are used in the analysis: 



• KQ85: The WMAP KQ85 cut with point source mask. Sky fraction used: 82% 

• KQ75: The WMAP KQ75 cut with point source mask. Sky fraction used: 72% 

• KQ75 ext.: The KQ75 cut extended with 5 degrees along the rim of the galaxy. Point 
source mask unchanged. Sky fraction used: 63% 

• |6| > 30 (sometimes referred to as the 60 degree cut): Same as the extended 
KQ75 mask, but with an additional 60 degree band cut on the galactic equator. Sky 
fraction used: 47%. 

• KQ85N equals the \b\ > 30 cut in the northern galactic hemisphere and the KQ85 cut 
in the southern galactic hemisphere. Sky fraction used: 65% 

• KQ85S equals the \b\ > 30 cut in the southern galactic hemisphere and the KQ85 cut 
in the northern galactic hemisphere. Sky fraction used: 64% 



3. Methodology: Power spectrum estimation on hemispheres 



3.1. Hemisphere spectra 



In ( lEriksen et al.l l2004al : iHansen et al.l l2004al ). we estimated the power spectrum on 
hemispheres centered on 164 di fferent positions on the sp here. In that analysis, we applied 
the Gabor transform approach (IHansen et al.ll2002l . |2003| ). In order to speed up the analysis 
without significant loss of precision, we will here apply the much faster MASTER algorithm 
(IHivon et al.ll2002l ). This allows the analysis of many more positions on the sphere as well 
as a large number of simulations for each WMAP year of observation and each channel. 



^http: / /lambda. gsfc.nasa.gov/ 
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The power spectrum is estimated from the pseudo power spectrum by (iHivon et al. 



20021) 



C(, = ^ K^g}{Ci' — Ni>) Ci 



t ~x ~Y 

2i+l 



where af^ is the spherical harmonic transform of channel X with a given mask, Ci is the 
corresponding pseudo power spectrum and Ne is the noise power sp ectrum. The coup ling 
kernel Ku/ depends on the mask applied to the data as detailed in ( IHivon et al.ll2002l ). In 
previous works we only used the auto-spectra, i.e. X = Y. Here we will also use the cross 
power spectra X 7^ F as an additional check. In the previous papers, we estimated the 
power spectrum on hemispheres centered on 164 different positions on the sphere. Here we 
use the positions of the 3072 pixel centers in a HEALPijfl A^sidc = 16 map. 

Using the above approach, we obtain for each multipole bin b a A^side = 16 map Mi{b) 
where the value of each pixel i corresponds to the i{£+ l)Ce power on a hemisphere centered 
on that pixel. In harmonic space, we bin the spectrum in 2 multipoles per bin such that 
if = 2b + 2 and ii = 26 + 3 are the first and last multipoles of bin b. We have also performed 
tests on more localized spectra, i.e. spectra estimated on disks of various sizes: 



Hemispheres (diameter 180°), 2 multipoles per bin, if = 2b + 2 and £1 = 2b + 3 
90° diameter disks, 4 multipoles per bin, if = Ab + 2 and £1 = Ab + 5 
45° diameter disks, 16 multipoles per bin, £f = 166 + 2 and £1 = 16b + 17 
22.5° diameter disks, 16 multipoles per bin, if = 166 + 2 and ii = 16b + 17 



In analogy to the tests considered in the previous papers ( Eriksen et al.ll2004al : lHansen et al 



2004al ). we tested the asymmetry for different multipole ranges constructing the following 
map: 



M 



Imirijlmax 



The max power spectrum ratio for a given multipole range was then defined as 



^lmin,lmax _ ^^^^ 



-M 



Iminjmax 



'M; 



Iminjmax ' ' 



where j is the pixel opposite to i. 



^http:/ /healpix.jpl. nasa.gov/ 
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Using the multipole ranges i = 2 — 19 and £ = 2 — 41 for which a significant asymmetry 
was found in previous papers, we still find a significant asymmetry with a similar direction 
and significance. 

The problem with this approach however, is that the max asymmetry axis as well as the 
significance is different for different multipole ranges. We made some attempts to determine 
whether the asymmetry continues to higher multipoles and if so, what is the highest multipole 
where the asymmetry is present. This turned out to be difficult because of the instability 
of significances and directions as a function of the multipole range. Whether we have a 
significant asymmetry or not depends on which multipole range we chose to look at. Also 
the axis of asymmetry is slightly different for different multipole ranges. 

We will now present an approach which, first of all will solve the problem of choosing 
which multipole ranges to look at, and secondly will tell us to which degree a complicated 
asymmetric model is preferred by the data rather than the isotropic model. The approach 



is inspired by a similar idea in iLand fc Magueijd (120071 ). Then in order to understand the 



results from the model selection method, we will present a second much simpler asymmetry 
test which is based on testing the alignment of the axes of asymmetry between independent 
multipole ranges. 



3.2. Model selection method 

We can look at the asymmetry in the following way: The asymmetry is the result of a 
dipole component in the maps which is common for several multipole bins b. The dipole 
for a given multipole bin b is given by 

^Im — 5^ ^Im ■ 

i 

We thus propose the following asymmetric model: 

4^ = a\^{0) + A{b)aim{9,(P) 

where a\^{0) is the random dipole expected in an isotropic model and aim{0,(f)) is the 
common dipole component. The parameters (0, 9) are assumed to be independent of the bin. 
We will allow variations of the asymmetry amplitude A with bin. To test this hypothesis, 
we will apply a simple fif? 

x^ = dtc-M, (1) 

where the elements of the data vector d are given by dmb = (Ami'^^^) ~ A{b)aim{d,(j)) for 
all m = [— 1, 1] and all bins b. Here 0^^(065) is the dipole of the map for the data to 
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be tested. The elements of the correlation matrix are simply Cmb,m'b' — {'Ami'^im')*) ■ 
will minimize this with respect to the parameters {A(b),6,(p). The resulting asymmetry 
direction is the best fit common dipole component for all bins included in the analysis. The 
value of the the minimum will be compared to the ior the isotropic hypothesis, 

Xl = dJC-Mo, (2) 

where the elements of do are given only by the observed dipole (io,m6 = ciimi'^^^)- The x^ 
improvement Ax^ = Xo~x'^ due to the additional parameters will be calibrated with a set of 
isotropic gaussian simulations. These simulations will give us the expected improvement Ax^ 
from the additional parameters to which the x^ improvement in the data can be compared. 
If the improvement Ax^ in the data turns out to be significantly better than in simulations, 
this would indicate that the asymmetric model is preferred by the data. 

We will test three different models for the asymmetry amplitude A{b), 

1. Constant amphtude: We will test (a) a three-parameter model (Aq, 9, 4>) with constant 
amplitude Aq for a given set of multipole ranges, (b) a four-parameter model where 
the asymmetry is assumed to have the constant amplitude Aq starting at £ — 2 until 
imax where the latter is the fourth free parameter, (c) a five-parameter model where 
both the lowest and highest multipoles in the asymmetric multipole range are free 
parameters. 

2. Linearly decreasing amplitude: In this four-parameter model we will assume the asym- 
metry to be maximal at i — 2 and decrease linearly according to A — Aq{1 — at) where 
a and Aq are free parameters. 

3. Gaussian multipole dependence: We will test a model where the asymmetry peaks at a 
certain multipole io and falls off to both sides following a Gaussian, A — Aoe^^^~^°^ 1^^'^ ^ . 
We will test (a) a four-parameter model which is assumed to peak at 1^ — 2 with a as 
free parameter and a five-parameter model where both £o and a are free parameters. 

Before presenting the results, we will show the procedure that we use to obtain these 
results. The following is the procedure that we use in the case of a 5 parameter model with 
a flat amplitude Aq (i.e. the amplitude is constant over a multipole range between fmin and 
^max and zero for all other multipoles). The free parameters in this model are thus: the 
two direction angles, the amplitude of the common dipole and the minimum and maximum 
multipoles of the multipole range where a common dipole component is found (i.e. where 
the amplitude is nonzero). The simulations are made with a maximum multipole of Lmax, so 
the maximum and minimum multipoles will be sought within the multipoles available from 
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the simulation in the range from L^\^ = 2 to Lmax- Since we will be looking for asymmetry 
extending over large multipole ranges, we will only look for multipole ranges with a minimum 
number of = ^max — ^min multipoles with a common dipole. 

1. We make two sets of 1400 WMAP simulations of a given band with a given mask. 

2. In each simulation, we estimate the power spectrum Cb{i) for a hemisphere centered 
on the center of each pixel i in the A^side = 16 HEALPix grid. We thus obtain 3072 
power spectra, one for each direction on the sky up to a multipole I = Lmax- For each 
simulations we thus have one map for each bin b. 

3. We extract the dipole from the map Mf for each bin b and each simulation. 

4. We use the first set of simulations to construct the correlation matrix C. 

5. We will now use the second set of simulations: For each simulation, we test all dif- 
ferent anisotropic models with different values for the five free parameters, imm, ^max, 
dipole amphtude Aq and direction {6, (/)). We will calculate the (equation [1]) of the 
simulated map for all these models and record the parameters of the model with the 
lowest x^- We use the following procedure to minimize the each simulation: 

(a) We start by making a loop over all possible multipole ranges between Lmin = 2 
and Ljnax which are the multipoles which we have available. If for instance -Lmax = 
300 and Ai = 200, we will go through all models with {imm = 2,£max = 201), 

(^min 2, £niax 221), (^niin '^j^ma.x 241), (^min 21,i'jnax 221), 

(^min 21, £jnax 241), (^min 101,i'max 300). 

(b) For each multipole range (^min, ^max), we make a loop over a large set of possible 
values for the model angle 9. 

(c) For each combination of imin, ^max and 6, we minimize the analytically with 
respect to the model parameters <p and ^o- Having found the values (p aiid 
that minimize the we now record the value of the minimum given these 
three values of the model parameters (^min, ^max, ^) as well as the best fit values 
Aq and (f) for this combination of (^min, ^max, ^)- 

(d) After terminating the loop over multipole ranges and directions 9, we now have 
a three dimensional array of minimum x^ values. We search for the minimum value 
of x^ in this array and thereby obtain the set of model parameters (^mim ^max, <; 
which minimizes the x^- We now compare this minimum x^ with the Xq (equa- 
tion [1]) obtained assuming the isotropic model. This difference Ax^ = Xo ~ is 
recorded for the given simulation. This value shows the improvement in for 
this given simulation when using a five-parameter anisotropic model. 
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6. After repeating the above procedure for all 1400 simulations, we repeat the same 
procedure on the WMAP data. We now have an array of improvements from a five- 
parameter model from all of the simulations as well as for the data. We now check the 
improvement of the data with respect to the simulations. If we quote a significance 
of 1% it means that 1% of the simulations had a similar or larger improvement of the 
using the five-parameter model. In the results we will also list the best fit multipole 
range (^min, ^max) as well as the amplitude A and direction {9, (j)) of the common dipole 
in the data. 

The exact procedure described above turned out to have convergence problems. The 
number of simulations used to obtain the correlation matrix was too small and caused small 
instabilities in the best fit parameters. 

The number of simulations necessary to get a converged covariance matrix would be too 
CPU demanding and we therefore chose the following solutions: By increasing the size of 
the power spectrum bins, the correlations between bins get smaller. Using power spectrum 
bins of A£ = 20, we found that the correlations can be ignored and only the diagonal part 
of the correlation matrix is included. With this bin size, identical results are obtained no 
matter whether the full correlation matrix or only the diagonal is used. With smaller bin 
sizes there was still some dependence on whether the full matrix was used or not. In the rest 
of this paper, we will therefore only use power spectra averaged in bins of 20 multipoles and 
a diagonal correlation matrix. Only in exceptional cases where we look at small multipole 
ranges will the bins of two multipoles be used. In this case, this will be clearly stated. 

In figure[T]we show the preferred direction for 1400 simulated maps using various galactic 
cuts. The plot shows the density of best-fit directions as a function of galactic co-latitude. 
The KQ85 and KQ75 cut show a uniform distribution of best-fit directions whereas for the 
larger cuts there is a preference for the poles. Thus we would expect that the preferred 
direction of asymmetry will be shifted away from the galactic plane for large sky cuts. 
Because of the large cut, power spectra estimated on hemispheres centered close to the poles 
will be similar in the polar area. This produces a dipolar structure with an axis pointing 
towards the poles. 



3.3. Test of alignment of multipole ranges 

In the results section we will see, using the model selection method described above, 
that a model with a common dipole component in the range from i = 2 — 600 is preferred 
by the data with high significance. There should thus be a strong correlation between 
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the distribution of power in different independent multipole ranges between £ = 2 — 600. 
Tlie purpose of the ahgnment test is to check if the asymmetry is distributed over the full 
multipole range and thus showing up as an alignment of the dipole of the power distribution 
of different independent multipole ranges. 

The idea for the alignment test is simple: we construct the maps 

b2 

M, (61,62) = Y,Mi{b), 

b=bl 

where we sum over (A) blocks of 20 multipoles and (B) blocks of 100 multipoles. Thus for 
test A, we obtain a set of maps Mi(2,21), Mi(22,41), Mi(42,61) etc. and for test B we 
obtain Mi{2, 101), Mi(102,201), Mi(202,301) etc. (due to different binning for 45 and 22.5 
degree tests, the ranges in test B will be as given in table [3]). 

For each map Mj(61, 62) the dipole is extracted (by a simple spherical harmonic trans- 
form on the map Mj(61, 62)) and the direction of the dipole is stored in a vector Vi where i is 
a multipole range (61, 62). In order to assess whether these directions for different multipole 
ranges are significantly more aligned in the WMAP data than in isotropic simulations, we 
define the mean angular distance 6 as 

9 = arccos Vt ■ Vj , 

where the sum over i and j is over subranges(61, 62) up to the maximum multipole for the 
given case. We will in the following quantify the alignment of the power distribution in the 
WMAP data by specifying the number of simulations with a lower mean angular distance 6 
between the dipoles of the power distribution. 



4. Results 

4.1. Results with model selection 

As a first test of the model selection approach, we used the constant amplitude model 
with £min and £max fixed at the ranges i = 2 — 19 and i = 2 — 41, allowing only the amplitude 
and direction of asymmetry to vary. In this case, we found the best fit asymmetry axis 
to be (138°, 220°) an d (108°, 227°) being within 10° - 15° of the asymmetry axis found in 
( Hansen et al. I l2004ah obtained with the method described in section [3m Note that due to the 



limited multipole range, we use power spectrum bins of two multipoles in this case, ignoring 
correlations between multipoles. In figure [2] we show the distribution of improvements 
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X^{isotropic) — ^^(Ao, (p, 6) with the three parameters (Aq, 6) obtained with 1400 isotropic 
simulated maps. The vertical line in the plots show Ax^ for the data. For £ = 2 — 41 we 
see that only 0.4% of the simulations have a drop in similar to the drop seen in the data, 
showing that the anisotropic model is actually preferred by the data. For the multipole range 
^ = 2 — 19 however, 30% of the simulations show a similar drop in and this asymmetry 
is therefore not significant alone taking into account the additional number of parameters 
required to describe it. 

We now allow first fmax (fixing ^min = 2) and then later also £min to be free parameters, 
i.e. we set the amplitude ^4 to a constant value Aq in a multipole range (£min,^max) and to 
zero for all other multipoles. Again we measured the drop in by the addition of 4 (fixing 
^min = 2 and varying £max) and 5 parameters (varying both ^min and ^max) for the simulations 
and compared to the data. 

In table [1] we show the results for larger scales. We extracted power spectra up to 
Lmax = 300 from WMAP simulation of different channels with different galactic cuts. As we 
are looking for asymmetric models extending over a large range in multipole space, we restrict 
our model search to models with at least 10 consecutive multipole bins (200 multipoles) with 
a common dipole component. The results for the 4-parameter fits (The four free parameters 
are the constant amplitude Aq, the direction {6, (j)) and the maximum multipole of asymmetry 
^max) are presented in the table. 

The table shows that within the multipole range i = 2 — 300, there is a dipole component 
with common amplitude and direction for the multipole range £ = 2 — 221. For the KQ85 cut, 
only 0.1 — 0.3% of the simulations show a similarly strong fit (similarly large improvement) 
for an asymmetric model. We see that the result is stable with frequency channel. The 
significance is dropping with larger galactic cuts, but even with the extended KQ75 cut, 
only 1.5% of the simulations give a similarly strong fit to an asymmetric model. Even with 
the largest cut, the direction is remarkably consistent, but note the expected shift away from 
the galactic plane as discussed above. The dipole fitting procedure has thus revealed that 
the hemispherical power asymmetry extends to at least £ = 221. 

The fact that the range i = 2 — 221 is found to be the best fit asymmetric range, 
does not mean that the asymmetry cannot extend beyond £ = 221. Our model consists 
of an isotropic field which in our analysis is considered noise and an anisotropic dipole 
component which is considered the signal. The isotropic 'noise' component to the dipole 
will randomly change the direction and amplitude of each single power spectrum bin away 
from the asymmetric direction. By including larger multipole ranges, this noise component 
is reduced and asymmetries extending over larger multipole ranges even beyond £ = 300 can 
give a good fit. 







50 100 150 

Co-latitude (degrees) 



Fig. 1. — The preferred direction for 1400 isotropic simulations using the KQ85 cut (sohd 
hne), KQ75 cut (dotted hne), extended KQ75 cut (dashed hne) and the 60 degree cut (dot- 
dashed hue). Here, the amphtude, direction and imax were free parameters and may thus be 
different for each simulation. 




Fig. 2. — Histogram of the improvements in for a three-parameter model with [6, (j), Aq) 
as free parameters for 1400 simulated maps using the WMAP V-band parameters with the 
KQ85 galactic cut. The left plot is for a model with asymmetry in the range £ = 2 — 19, the 
right plot is for a model with asymmetry in the range ^ = 2 — 41. The horizontal line is the 
improvement in the foi" the 5 year WMAP data. 
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From table [T] we also see that the best fit asymmetric multipole range is slightly larger 
for the larger galactic cuts. For KQ85 £ = 2 — 221 is the best fit asymmetric range whereas 
for KQ75 i = 2 — 281 gives a better fit and for larger cuts even i = 2 — 300 is preferred. As 
discussed above, the fact that a smaller range i = 2 — 221 gives the best fit for KQ85 does 
not mean that the larger range i = 2 — 281 is a bad fit for this mask. In the table we have 
included in parenthesis the significance for £ = 2 — 281 for KQ85. Clearly i = 2 — 281 is also 
significantly asymmetric for KQ85, but because of the random noise component, £ = 2 — 221 
gives a shghtly better fit. 

Note that we have not included the results of the 5-parameter fits. The reason for this 
is that the 5-parameter fit in all these cases show identical best-fit model parameters (that 
is, they show that the best fit value for the fifth parameter ij^^^ equals £min = 2 for all cases) 
to the 4-parameter fits. The significances for the 5-parameter model are generally lower 
than the 4-parameter model even if the best fit parameters are the same. The reason for 
this is that significances for the 4-parameter model is calibrated with improvements in 
4-parameter fits to isotropic simulations whereas the 5-parameter model is calibrated with 
improvements in 5-parameter fits to simulations. The simulations will in general have 
a larger improvement for a 5-parameter model than a 4-parameter model. The data 
however will have the same for a model with exactly the same parameters. There will 
therefore be more simulations with a larger improvement for the 5-parameter fit than for 
the 4-parameter fit even if the of the data is the same. 

To investigate properly the maximum multipole for asymmetry, we run a set of simula- 
tions with I/max = 500 and Lmax = 800. In order to reduce the CPU time for the Lmax = 800 
case, we now use disks of 90, 45 and 22.5 degree diameter. In this way we also obtain more 
localized spectra. For the 45 and 22.5 degree disks, we needed to use multipole bins of 16 
multipoles in the power spectrum estimation, instead of 20 used in the previous analysis. 
For the 22.5 degree disks, the variance of the power spectrum estimate close to the galactic 
plane was so large that the map needed to be normalized by its standard deviation 
(obtained from simulations) before a dipole fit could be performed. As a result, the values 
for the amplitude obtained in this case is different from the amplitudes obtained for other 
disk sizes. 

The results are shown in table [2l The hemisphere results with Lmax = 500 show that the 
asymmetry extends to £ = 481 with a similar direction of asymmetry as for lower multipoles 
but now with a lower significance {p = 2.3%). For the more localized spectra however, the 
asymmetry is found highly significant (p = 0.4% for the most localized spectra) for the range 
i = 2 — 600 for all disk sizes, but no evidence is found for an asymmetry extending beyond 
i = 600. In the same table we also show results for the dipole fit on smaller disks including 
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only multipoles up to Lmax = 300 in order to check for consistency with the hemisphere 
results. Clearly a consistent dipole fit is found also in this case with more localized power 
spectra. 

Finally we will make a consistency check by performing the 5-parameter dipole fit in 
individual subranges of 100 multipoles from i = 2 to i = 800 using 22.5 degree disks. Note 
that the size of these subranges is not exactly 100 because each power spectrum bin has 16 
multipoles. In table [3] we show the results. We see clearly that the best fit direction in each 
subrange up to £ = 600 is consistent with the best fit direction {6 = 107, = 216) for the 
full range £ = 2 — 591 for this disk size. The two bins above £ = 600 however show a very 
different dipole direction. We see that the asymmetry can be seen as an alignment of the 
power distribution dipoles between multipoles from £ = 2 up to ^ = 600. We will now study 
this alignment in more detail. 



4.2. Results with the alignment test 

Before looking at significances, we will illustrate the direction of the dipoles of individual 
100 multipole blocks with some figures. In figure [3] we show the distribution of power in the 
WMAP V+W band data using hemispheres (KQ85 cut was used in the power spectrum 
estimation). Each map shows the distribution of power Mj(61,62) for a given 100-multipole 
range. We see already by eye that there is a clear dipolar distribution and that the direction 
of the dipole is very similar in each case. In figure H] we show the position of the dipole for 
subranges of 100 multipoles. The color of the disk indicates the multipole range (see table [3] 
for the exact ranges used). The results in this plot is taken from power spectrum estimates on 
disks with diameter 22.5 degree using the KQ85 galactic cut. We see that all the individual 
multipole ranges have dipoles pointing in a direction close to the best fit dipole for the full 
range £ = 2 — 600 indicated by the white hexagon. 



In [Hansen et al.l (l2004a] ). we point out that the outliers in the full sky spectrum at 
£ = 22 and ^ = 40 seem to be associated with the asymmetry: The high outlier at ^ = 40 
was associated with the high power in the hemisphere of maximum asymmetry and the low 
outlier at £ = 22 was associated with the low power in the opposite hemisphere and. In figure 
|5] we show the distribution of power in these two bins as well as for the first bin £ = 2 — 3 
and the bin i = 28 — 29 which is a particularly asymmetric bin. The conclusions from 



Hansen et al.l (l2004a( ) still hold. One can see that the same dipolar distribution of power 



which is also seen in the 6 subranges between £ = 2 — 600 in figure O 



In table H] we show the results for the alignment test using blocks of 20 multipoles. 
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Fig. 3. — The distribution of power in blocks of 100 multipoles estimated on hemispheres 
for the combined V+W band using the KQ85 sky cut. Note the similarity with the single 
multipole bins in figure [5] 
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Table 1. The table shows significances (in %) and parameters (fmaxi (f), Aq) of the best 
fit asymmetric 4-parameter model. Approximate Fisher matrix error for 6, (f) and Aq are 

also given. Please refer to the text for details about the asymmetric models and their 
parameters. The significances specify the percentage of simulated maps with a larger drop 
in for the asymmetric model (considering only models with A£ > 200) than found in the 
WMAP data. The results are based on 1400 simulations. 



Channel 


mask 


^-range 


•''max 


^(deg) 


0(deg) 


p(%) 


AO 


A(f) 


Ao xlO"^ 


Q 


Kq85 


[2,300] 


221 


104 


226 


0.1 


10 


10 


1.6 ± 0.4 


V 


Kq85 


[2,300] 


221(281) 


107 


226 


0.3(0.3)* 


10 


10 


1.5 ± 0.4 


w 


Kq85 


[2,300] 


221 


103 


229 


0.1 


9 


9 


1.7 ± 0.4 


V 


Kq75 


[2, 300] 


281 


112 


216 


3.1 


14 


13 


1.2 ± 0.4 


V 


Kq75 ext. 


[2, 300] 


300 


114 


202 


1.5 


12 


13 


1.3±0.4 


V 


\b\ > 30 


[2,300] 


300 


131 


170 


17 


18 


19 


1.0 ± 0.4 



* This results was obtained with 4200 simulations 




Fig. 4. — The directions of the dipoles of the power distribution in blocks of 100 multipoles 
estimated on disks with diameter of 22.5 degree for the combined V+W band using the 
KQ85 sky cut. The colour of the disks indicate the center of the given multipole range. The 
white hexagon indicates the best fit dipole direction for the full range i = 2 — 600. The 
ecliptic poles are indicated by crosses. 



-17- 



Table 2. The table shows significances (in %) and parameters (£max, (f), Aq) of the best 
fit asymmetric 4-parameter model. Approximate Fisher matrix error for 6, and Aq are 

also given. Please refer to the text for details about the asymmetric models and their 
parameters. The significances specify the percentage of simulated maps with a larger drop 
in for the asymmetric model than found in the WMAP data. The results are based on 
1400 simulations with the co-added V+W channels. 



Diameter 


mask 


£-range 




^(deg) 


0(deg) 


p(%) 




A0 


Ao xlO-^ 


Ai > 400 




180° 


Kq85 


[2, 500] 


481 


102 


235 


2.3 


12 


12 


0.7 ± 0.2 


180° 


Kq75 


[2, 500] 


481 


104 


224 


14 


16 


16 


0.5 ± 0.2 


90° 


Kq85 


[2, 800] 


601 


105 


225 


2.6 


11 


11 


1.3±0.3 


45° 


Kq85 


[2, 800] 


591 


102 


223 


0.6 


9 


10 


1.7±0.4 


22.5° 


Kq85 


[2, 800] 


591 


107 


216 


0.4 


11 


10 


1.4±0.4 


AC > 200 




90° 


Kq85 


[2, 300] 


221 


102 


229 


0.3 


11 


11 


2.6±0.7 


45° 


Kq85 


[2, 300] 


223 


98 


227 


1.1 


11 


12 


3.0±0.9 


22.5° 


Kq85 


[2, 300] 


223 


97 


220 


2.1 


14 


14 


0.21±0.07 




Fig. 5. — The distribution of power i{i + l)Ce for the multipole ranges i = 2 — 3 (upper 
left), £ = 22 - 23 (upper right), £ = 40-41 (lower left) and £ = 28 - 29 (lower right) 
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We limit the maximum multipole to 300 because noise is getting important after ^ = 300 
increasing the variance of the directions of blocks with just 20 multipoles. Later for blocks 
of 100 multipoles we will also consider higher multipoles as the noise is reduced in each 
block when averaging over 100 instead of 20 multipoles. In the table we have considered 
the alignment test for the V band using different galactic masks. The numbers indicate 
the percentage of simulations with a lower mean angle 9 between dipole directions. For 
£ = 2 — 300 there is a significant (~ 1% level) alignment between dipole directions for the 20 
multipole blocks. Note in particular that the for the extended KQ75 cut, none of the 1400 
simulations are as strongly aligned as the WMAP data. Thus, the asymmetry is strong even 
with a large galactic cut. We also see that there is a significant (2 — 4% level) alignment when 
considering only the first five 20-multipole blocks on large scales £ = 2 — 100. Considering 
only the blocks on small scales i = 200 — 300 excluding the first 100 multipoles however, 
there is no significant alignment present. 

In table Owe show the results from the alignment test of 100 multipole blocks using more 
localized spectra. The V+W map with the KQ85 galactic cut was used in this analysis. We 
clearly see the result of the strong alignment which was already obvious in table [31 Using 
more localized spectra the alignment appears even stronger than for hemisphere spectra. 
Note that for the range £ = 2 — 600 none of the simulations show a similarly strong alignment 
for any of the disk sizes. In particular for the 22.5° (diameter) disk results, none of 9800 
simulations have a similarly strong alignment which is close to a 4a detection of asymmetry. 
Note also that the alignment is highly significant also for separate multipole ranges at small 
and large scales, for instance £ = 2 — 300 and £ = 300 — 600. 

In figure [6] we show the spectra in the best fit dipole direction for £ = 2 — 600 for various 
disk sizes from hemispheres to 45° (diameter) disks. Also in these plots we see that the 
difference between the spectra in the opposite directions becomes larger with more localized 
spectra. In particular the first part of the spectrum £ = 2 — 100 as well as the amplitudes of 
the first two peaks are clearly different. We have investigated whether also the positions and 
not only the amplitudes of the first two peaks may have a similarly asymmetric distribution 
on the sky. For the localized spectra obtained on 90° (diameter) disks using the V+W map 
with the KQ85 cut, we made a fit to the first two peaks. In figure [7] we show the distribution 
of the multipole position of the first and second peak. No dipole structure similar to what 
was found for the power spectrum amplitude is seen in these figures. 
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Table 3. The table shows significances (in %) and parameters (^mm, ^max, ^5 (f'j ^o) of the 
best fit asymmetric 5-parameter model. Approximate Fisher matrix error for 9, <p and Aq 
are also given. Please refer to the text for details about the asymmetric models and their 
parameters. The significances specify the percentage of simulated maps with a larger drop 
in for the asymmetric model than found in the WMAP data. The results are based on 
1400 simulations for the co-added V+W channels. 



Diameter 


mask 


£-range 


^min 




^(deg) 


0(deg) 


p(%) 




A0 


Ao xlO-^ 


Ae > 40 




22.5° 


Kq85 


[2,95] 


2 


63 


110 


226 


8.4 


15 


16 


0.38±0.14 


22.5° 


Kq85 


[96,191] 


96 


191 


100 


200 


19 


16 


16 


0.27±0.10 


22.5° 


kq85 


[192,303] 


208 


281 


100 


238 


99 


35 


34 


0.13±0.11 


22.5° 


kq85 


[304,399] 


352 


399 


83 


182 


77 


24 


24 


0.23±0.13 


22.5° 


kq85 


[400,495] 


432 


479 


113 


224 


36 


19 


19 


0.32±0.14 


22.5° 


kq85 


[496,591] 


496 


591 


112 


210 


48 


20 


20 


0.21±0.10 


22.5° 


kq85 


[592,703] 


608 


687 


45 


111 


43 


20 


24 


0.26±0.11 


22.5° 


kq85 


[704,799] 


736 


799 


35 


47 


63 


26 


36 


0.25±0.13 



Table 4. For each entry in the table we have calculated the mean angle 9 between the 
dipole directions for all blocks of 20 multipoles within the given multipole range (refer to 
the text for details of how 9 is calculated). The numbers given in the table is the percentage 
of simulated maps with a lower mean angle 9. The results are based on 1400 simulations. 
Zero entries means that none of the simulated maps had a similarly low mean angle. 



mask: 


KQ85(V) 


KQ75(V) 


KQ75ext(V) 


|6| > 30(V) 






£ = 2 - 300 


0.9 


0.9 





5.6 


£ = 2 - 200 


1.8 


5.5 


0.8 


13 


£ = 2 - 100 


2.0 


4.3 


3.4 


16 


£ = 100 - 300 


9.0 


7.4 


1.8 


35.5 


£ = 200 - 300 


67.5 


24.6 


20 


49.9 
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Table 5. For each entry in the table we have calculated the mean angle 9 between the 
dipole directions for all blocks of 100 multipoles within the given multipole range (refer to 

the text for details of how 6 is calculated). The numbers given in the tabic is the 
percentage of simulated maps with a lower mean angle 6. The results are based on 1400 
simulations except for the results for 22.5° disks which is based on 9800 simulations. Zero 

entries means that none of the simulated maps had a similarly low mean angle. The 
combined V+W map was used in obtaining all results in this table. The KQ85 mask was 

used when other mask is not specified. 



disk size: 


180° 


90° 


45° 


22.5° 


90°(KQ85N) 






£ = 2 - 800 




6.1 


8.3 


9.2 


53 


£ = 2 - 700 




0.5 


0.5 


0.4 


15 


£ = 2 - 600 













1.1 


£ = 2 - 500 


0.1 








0.04 


3.9 


£ = 2 - 400 


0.1 





0.1 


0.3 


4.6 


£ = 2 - 300 


0.5 


0.6 


0.7 


1.8 


11 


£ = 2 - 200 


14 


13 


9.2 


18 


41 


£ = 200 - 600 













2.7 


£ = 300 - 600 




0.1 


0.1 


0.1 


7.3 


£ = 400 - 600 




3.6 


0.4 


0.5 


6.6 
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Fig. 6. — Power spectra obtained in opposite directions. The upper two plots are spectra 
estimated on hemispheres, one centered at (102°, 235°) (black line) and one in the opposite 
direction (grey line). The left plot shows the spectra after foreground subtraction (by the 
WMAP team), the right plot shows the same spectra before foreground subtraction. We 
see that if the asymmetry is larger than the foreground correction which is applied to the 
maps. The lower two plots show the more localized spectra taken in opposite directions, left 
plot for spectra estimated on 90 degree (diameter) disks, right plot for 45 degree disks. The 
asymmetry is more pronounced for more localized spectra. 
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Position of first peak 




3.1e + 03 ^^^^^^^^^^K^ ^^^^^^^H <:.3e + 03 central multipole 



Position of first peali 




Fig. 7. — For each position on the sphere, the colour indicates the multipole positions of 
the first (upper plot) and the second (lower plot) Doppler peak. The spectra were estimated 
on disks with diameter of 90° centered on the given position on the sphere. The combined 
V+W map with the KQ85 sky cut was used for the power spectrum estimation. 
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4.3. Testing foregrounds and systematics 

In this section we will perform several tests in order to investigate whether foreground 
residuals or instrumental systematic effects may cause the observed asymmetry. In particu- 
lar, we will look at the cross spectra instead of the auto-spectra, we will look at the WMAP 
data year by year and finally we will study in detail whether there are still clear signs of the 
asymmetry outside the \b\ > 30 cut. 

The cross spectra based on aim obtained from different channels and years of obser- 
vations are_Jess_groiieto_systematical errors and in particular to uncertainties in the noise 



model ( iHinshaw et al.ll2003l ). We have obtain the hemisphere spectra based on the spectrum 
obtained as a mean of all 780 possible combinations of the channels Q, V and W as well as 
the five years of observation. Making simulations with all 780 cross-spectra turned out to 
require too much CPU time and we were therefore not able to perform a full statistical test 
using the cross-spectra. Using the WMAP data alone we found that the direction of the 
dipole using the maps based on cross-spectra is consistent with the dipole based on the 
auto-spectrum. 

Similarly we have studied the direction of the dipole for each single year of observation. 
We found that these directions are consistent with the direction obtained with the co-added 
maps with all years included. There is thus no sign of systematic errors in specific years 
causing the asymmetry. We have also considered the difference between the power distri- 
bution maps Mi{bl,b2) obtained with different channels. Foreground residuals causing the 
asymmetry would show up in these differences between channels. The dipole directions of 
the two difference maps Q-V and V-W are not similar to the direction of asymmetry. There 
is thus no sign of a frequency dependent foreground with the dipolar power distribution 
which we have detected in the individual bands. 

As discussed above, the asymmetry is no longer significant for the \b\ > 30 cut, but 
the fact that the direction of asymmetry is still (within the error bars) consistent with 
the direction of asymmetry found for smaller galactic cuts is a strong argument against 
foreground residuals causing the asymmetry. 

Still, in order to make sure that the drop in significance for the 60 degree cut is not due 
to the fact that the large mask is excluding some galactic residual causing the asymmetry, 
we made some further tests. First we made a mask which was equal to the regular KQ85 
cut for the southern galactic hemisphere, but had an extended 30 degree galactic cut in 
the northern galactic hemisphere (this is the KQ85N mask). We also made a similar mask 
extending only in the southern hemisphere (this is the KQ85S mask). 

Using the southern KQ85S mask (for i = 2 — 300), the significance of the asymmetry is 
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still high [p = 1.2%) for the range i = 2 — 261 (the direction of asymmetry is consistent with 
the results above). For the northern KQ85N mask the significance has dropped to 12%, but 
the direction of asymmetry is still consistent. Using the more localized spectra estimated on 
90° disks, we see from table that the alignment between 100-multipole blocks is significant 
at the 1.1% level for £ = 2 - 600 using the KQ85N cut. For ^ = 2 - 400 the alignment 
is still significant at the 2a level. This result combined with the fact that the direction of 
asymmetry has changed little with the large \b\ > 30 cut as well as the consistency of results 
using different frequency bands shows that an explanation of the asymmetry in terms of 
foreground residuals is difficult to make consistent with the results presented in this paper. 



4.4. Non-fiat amplitude 

We have so far discussed a model where the common dipole has the same amplitude 
Aq over the full multipole range in the best fit model. We have also tested a model with A 
decreasing linearly with multipole as well as a Gaussian shape of the amplitude. Both these 
models are described in detail in section [3^ In table [6] we present the results. The first part 
of the table shows the result with the linear fit where a is the parameter describing how fast 
A decreases with multipole. We see that a model with decreasing amplitude is preferred by 
the data. In figure [H] we show the multipole dependence of the asymmetry. The asymmetry 
decreases and vanishes close to £ = 600 consistent with the results in the previous chapters. 

In the lower part of the table, the results with the Gaussian model are shown. We shows 
the results for the 4-parameter model where the peak of the Gaussian is forced to ipeak = 2 
with the width a allowed to vary as well as for the 5-parameter model with £peak as an 
additional free parameter. With the exception of the 45 degree disk results, the 5 parameter 
model finds the same best fit model with ipeak = 2 as the 4-parameter model. In figure [8] we 
have plotted the best fit Gaussian model for 90 degree disks on top of the linear model. We 
see that the two models show a consistent decrease in the amplitude of the asymmetry. We 
conclude that the asymmetry is larger for smaller multipoles and decreasing continuously 
towards i = 600 where it disappears. 



5. Conclusions 



We have re assessed the asymmetr y in t he distribution of CM B fluctuation power on the 
sky reported in lEriksen et al.l ( j2004al ) and iHansen et al.l (j2004al ) . In order to test whether 
an anisotropic model of the CMB fluctuations is actually preferred over an isotropic model 
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Table 6. Significances and directions of tlie common dipole component for asymmetric 
models with linearly decreasing or Gaussian amplitude profile. We show the significances 
(in %) and parameters {6, (j)) as well as a for the linear model. For the Gaussian model we 
show the best fit parameters a (included only in the 5-parameter fit) and ipeak- Please refer 
to the text for details about the asymmetric models and their parameters. The significances 
specify the percentage of simulated maps with a larger drop in for the asymmetric 
model than found in the WMAP data. The results are based on 1400 simulations. 



Channel 


Mask 




9 (deg) 


(f) (deg) 


a 


p(%) 


Linear 




V+W 


KQ85 


0.0018 


99 


229 




0.6 


v+w 


KQ75 


0.0018 


100 


218 




17 


V+W 


KQ85 (90° disks) 


0.0015 


100 


230 




1.3 


v+w 


KQ85 (45° disks) 


0.0014 


99 


228 




1.4 


Gauss 




V+W 


KQ85 


NA 


230 


100 


228 


0.4 


v+w 


KQ85 


2 


230 


100 


228 


1.7 


v+w 


KQ85 (90° disks) 


NA 


252 


100 


230 


0.8 


v+w 


KQ85 (90° disks) 


2 


252 


100 


230 


14 


v+w 


KQ85 (45° disks) 


NA 


312 


100 


228 


1.4 


v+w 


KQ85 (45° disks) 


408 


312 


100 


228 


9.1 




0.0 L , , I , , , I , , , ''V-.- 

200 400 600 800 

Multipole 

Fig. 8. — The multipole dependence of the amphtude of asymmetry A. We show the results 
of the fit to a linear and Gaussian model of A. The results were obtained from spectra 
estimated on 90 degree (diameter) disks using the V+W band with the KQ85 sky cut. 
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taking into account the additional parameters required, we implement a new model selection 
procedure. We model the asymmetric distribution of power on the sky as a dipole in the 
power distribution. Note that this is not to say that the CMB fluctuation field has a dipole, 
rather that the power for a certain scale (multipole) has a dipole distribution on the sky. We 
use a model where there is a common dipole component in the power distribution for a set 
of multipoles in a range [^min,^max] where ^min and £max are free parameters as well as the 
direction (6*, </>) and the amplitude Aq of the dipole. We use a approach to find the best 
fit model parameters among these 5 parameters. 

We first investigated a model where we assume the asymmetry to start at £mm = 2, 
reducing the number of free parameters in the model to 4. Using power spectra estimated 
on hemispheres in the combined V+W band with the KQ85 galactic cut, we find a strong 
asymmetry in the multipole range £ = 2 — 221 with an axis pointing in the direction {9 = 
107° ± 10°, (I) = 226° ± 10°) (which is the direction where the power is largest). Only 0.3% 
of the simulated isotropic maps show a similarly strong asymmetry. Performing the same 
test on the Q, V and W bands individually as well as with different galactic cuts, a similar 
asymmetry is found. The significance is reduced for larger cuts but for an extended KQ75 
cut (excluding 37% of the sky), only 1.5% of the simulated isotropic map show a simiarly 
strong common dipole component in the range £ = 2 — 300. Using more localized power 
spectra estimated on smaller disks, we perform the same tests on the V+W band using the 
KQ85 cut. Smaller disks allow faster spherical harmonic transforms and allows the analysis 
to include multipoles up to £ = 800. We find that the range i = 2 — 600 is asymmetric with 
dipole direction {9 = 107° ± 11°, = 216° ± 10°) using the smallest 22.5° diameter disks. 
Only 0.4% of the simulated maps show similar asymmetry. In figure [6] we show the spectra 
in the two opposite parts of the sky. The spectra are clearly different for the largest scales 
as well as around the two first peaks. 

Including ^min as a free parameter, the same anisotropic model is favored and the best 
fit value for the first multipole with a common dipole is £min = 2. We therefore concluded 
that the 4 parameter model was sufficient to describe the asymmetry. Testing models with 
a multipole dependent dipole amplitude, we found that a model where the asymmetry is 
maximum for small £ and decreasing with increasing £ vanishing at about £ ~ 650 gives a 
good fit to the data (0.4%). 

To check whether the asymmetry is present in the full range £ = 2 — 600 or only for some 
multipoles, we performed a second simpler test of asymmetry. We found the dipole direction 
for the power distribution in multipole ranges of 100 multipoles, £ = 2 — 101, £ = 102 — 201, 
etc. to £ = 502 — 601. We thus obtained 6 dipole directions from 6 independent multipole 
ranges. The power distribution in these ranges are shown in figure O We found that these 
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6 dipoles were much more aligned in the WMAP data than in isotropic simulations. In 
fact, using power spectra estimated on 22.5° diameter disks from the combined V+W bands 
with KQ85 galactic cut we found that none of our 9800 simulations show a similarly strong 
alignment between these 6 dipoles. The dipole directions for these 6 ranges are shown in 
figure m We find that the spatial distribution of CMB fluctuations is strongly correlated 
between small and large angular scales. 

The fact that all the frequency channels, all years of observations and also tests using 
cross power spectra show a similarly asymmetric distribution strongly disfavors an expla- 
nation in terms of systematic effects and residual galactic foregrounds. Further, in the 5 
year WMAP data a different approach to foreground subtraction than for the first year data 
was applied, still the asym metry remained sign ificant at the same level for the large scale 



asymmetry £ = 2 — 40. In (IHansen et al.l |2006| ) we also showed that a blind approach to 



foreground subtraction did not change this result for large scales. 

Our results indicate that the reported common asymmetric axis extending over a large 
range in scales is highly unlikely to be a statistical fluke. Foregrounds and systematic effects 
do not seem to be probable explanations. The CMB does seem to have an uneven power 
distribution on the sky over a large range of angular scales. An important task for further 
research is to find a physical explanation for this asymmetry which can predict possible 
effects on CMB polarization to be tested in future experiments. 
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